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Two-Level  Supply  Chain 


•  The  buyer  randomly  receives  demand  from  its  customers  and 
places  orders  of  size  Q  from  the  supplier  when  its  inventory  level 
drops  below  reorder  point  R. 


Inventory 


' 

Demand 

(Random) 

Supplier 

Order 

Buyer 

•  The  supplier  receives  these  orders  and  ships  inventory  to  the 
buyer.  The  supplier  orders  a  quantity  from  its  supplier  in  integer 
multiples  N  of  Q. 
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Two-Level  Supply  Chain 


•  Buyer  cost  function: 

TCb(Q,  R,  V)  =  (Kb-V)~  +  hb-(0.5Q  +  R-E<iX))  +  7r-SR(R)- 

i — >•  Terms  are  annual  ordering,  holding,  and  stockout  costs;  V 
is  a  per  order  rebate  coordination  incentive  (Cobb  and  Johnson 
2013). 


Expected  shortage  per  cycle:  Sr(R )  = 


R )  •  fx(x )  dx 


>  fx  is  the  probability  density  function  (PDF)  for  lead  time 
demand. 
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<0  ^ 


Two-Level  Supply  Chain 


•  Supplier  cost  function: 


TCS(Q,  N,  V ) 


+  hs(N  -  1)0.5Q 


i— >  N  =  integer  multiple  of  buyer’s  order  size  Q 
i — *  +  =  annual  ordering  costs 

i — >  hs(N  —  1)0. 5Q  =  annual  holding  costs 
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Lead  Time  Demand 


•  LTD  follows  a  compound  probability  distribution.  Suppose  L 
is  a  random  variable  for  LT  and  D  represents  random  DPUT. 
LTD  is  a  random  variable  X  determined  as 

X  =  D\  +  L>2  +  Di  Dl  . 

Therefore,  X  is  a  sum  of  random,  independent  and  identically 
distributed  (i.i.d.)  instances  of  demand.  The  mean  (nx)  ancl 
variance  (crj)  of  X  can  be  calculated  as 

E(X)  =  E(L)-E(D )  and  Var(X)  =  E(L)-Var(D)  +  [E(D)]2 -Var(L)  . 
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Example  (Eppen  &l  Martin  (1988)) 

•  Daily  demand  is  normally  distributed:  Di  ~  N{ 40,30). 

•  Lead  time  (in  days)  is  discrete:  Q.L  =  {7,12,14,15,16,25} 
each  with  probability  1/6  ( E(L )  =  14.8;  Var(L)  —  29.1). 

•  In  this  case,  E(X )  =  fix  =  14.8  •  40  ss  593  and 

Var(X)  =  (j\  =  E(L)  ■  Var(D)  +  [E(D)]2  ■  Var(L)  ss  47000 


•  All  previous  methods  for  setting  and  N*  assume  LTD 

is  normal;  Eppen  and  Martin  (1988)  demonstrate  calculation  of 
a  service  level  in  a  single-firm  context. 
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Example  (Eppen  &l  Martin  (1988)) 

•  The  actual  LTD  distribution,  fx,  is  a  mixture  of  normal  distri¬ 
butions. 
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Two-Level  Supply  Chain 


•  Research  issues  to  be  addressed: 

1.  Analytical  solutions  for  Q*t  R*,  and  N*  in  the  two-level  supply 
chain  problem  assume  (and  require)  normality. 


2.  Methods  for  modeling  LTD  distributions  often  make  unreal¬ 
istic  distributional  assumptions. 


3.  In  practice,  the  actual  LT  and  DPUT  distributions  are  likely 
unknown  -  the  solution  here  uses  empirical  data. 


i — >■  For  the  example,  suppose  a  modest  amount  of  historical  data 
is  available  on  daily  demand  and  lead  time  on  previous  orders. 
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Example  (Eppen  &i  Martin  (1988)) 

i — ^  500  observations  of  daily  demand  with  d  =  39.66  and  5 
30.64. 


•  This  is  a  random  sample  from  the  7V(40,30)  distribution. 
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Example  (Eppen  &  Martin  (1988)) 

i — ^  50  observations  of  lead  time  with  I  =  10.8  and  sj  —  12.52. 


•  This  is  a  random  sample  from  the  discrete  LT  distribution. 
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Constructing  the  LTD  Distribution 

i — ^  500  observations  of  daily  demand  with  d  =  39.66  and  5 
30.64 
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•  These  are  assumed  to  be  i.i.d.  observations. 
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Constructing  the  LTD  Distribution 


The  most  likely  empirical  LT  value  is  7  days 


7-day  total 
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Sum  daily  demand  over  each  7-day  period  in  the  dataset 


Constructing  the  LTD  Distribution 

i — ^  The  lead  time  demand  dataset  given  LT  of  7  days  with  Nj  = 
71 
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•  Fit  a  mixture  of  polynomials  (MOP)  distribution  (Shenoy  2012) 
to  this  data.  This  distribution  will  be  the  approximate  LTD 
distribution  conditional  on  L  —  7,  or  fx\L=7- 

•  Similar  distributions,  fx\L=£>  wil1  iDe  constructed  for  t  —  12, 14, 
15, 16,25. 


13 


Mixture  of  Polynomials  (MOP) 


i — *  The  approximate  LTD  distribution  conditional  on  L  =  7,  or 

fx\L—7 


fx\L=l(.x)  ~ 

I  -3.3726  +  0.0241a;  -  0.000043a;2  236  <  a;  <  271 
I  0.7102  -  0.0061a;  +  0.000013a;2  271  <  x  <  306  . 


•  The  MOP  has  n  —  2  pieces  and  is  degree  d  —  2  (or  is  third 
order). 

i — *  The  MOP  was  constructed  using  a  linear  combination  of  13- 
spline  functions  (Lopez-Cruz  et  al.  2012). 
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Mixture  of  Polynomials  (MOP) 


i — >■  fx\L=7  is  a  mix'ture  of  four  B-spline  functions 


•  The  B-splines  are  defined  recursively  based  on  the  split  points 
in  the  domain.  Mixing  coefficients  are  determined  via  maximum 
likelihood  (Zong  2006). 
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Mixture  of  Polynomials  (MOP) 

i — >  The  approximate  LTD  distribution  conditional  on  L  =  7,  or 
fx\L=7  overlaid  on  the  N( 7  •  40, 7  •  30)  distribution. 


•  Recall:  the  MOP  is  not  fit  to  the  normal  PDF,  but  rather  a 
small  sample  of  data  generated  from  the  normal  PDF. 
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Mixture  of  Polynomials  (MOP) 


•  Selection  of  d  and  n  for  B-spline  estimation  is  a  trade-off  be¬ 
tween  accuracy  and  complexity;  higher  values  can  also  lead  to 
over-fitting. 

i — >  Select  d  and  n  to  maximize  the  Bayesian  information  criterion 
(BIC): 

BIC  (fx\L=e  O)  ,V)=C  (v\fx]L=e  (x))  -  ((m  -  1)  log  N)/2  . 

The  second  term  is  a  penalty  for  adding  parameters  to  the  model 
(£  is  the  likelihood  of  the  data  given  the  model). 

i — ^  In  practice,  once  we  settle  on  d  and  n,  we  may  not  go  through 
this  step  (n  —  2  &l  d  —  3  maximized  BIC  in  this  example). 
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Mixture  of  Polynomials  (MOP) 


•  What  does  “over-fitting"  look  like?  Left:  n  —  28^.d  —  3]  Right: 
n  —  3  &l  d  =  5. 


•  Left:  c(v\fL=7)  =  -284.0  &  BIC:  -290.4  ;  Right:  c(v\ fL=1) 
—  —280.9  BIC:  —293.7  (Bonus:  model  on  left  entails  less 
computational  complexity). 
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Mixture  of  Polynomials  (MOP) 

i — *  Fitting  process  repeated  for  t  —  12,14,15,16,25. 


The  approximate  LTD  distribution: 

fx 0)=  E  P(L  =  £(k))-f: 

fO)=i 


1000 


L=£(k)(X) 


Mixture  of  Polynomials  (MOP) 


i — >  fx  is  overlaid  on  the  mixture  of  normal  distribution  fx  (the 
“actual"  distribution),  but  fx  was  not  generated  using  knowl¬ 
edge  of  the  underlying  distribution. 
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Mixture  of  Polynomials  (MOP) 


i — *  fx  is  relatively  compact  -  it  has  15  pieces  and  is  a  2nd  degree 
polynomial. 
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Mixture  of  Polynomials  (MOP) 


i — >  How  is  this  useful? 

Recall:  Expected  shortage  per  cycle:  Sr(R )  = 


*x 


max 


'R 


(x— R)-fx(x)  dx 


i — ►  SR  can  be  calculated  in  closed-form  (here  a  10-piece,  5th 
degree  MOP),  so  TCb  is  closed-form. 
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Finding  Optimal  Policies 

Test  case  (from  CH):  Kb  =  50,  Ks  =  150,  hb  =  5,  hs  =  12.5, 
7r  =  6,  250  working  days  so  Y  —  250  *  40  =  10000. 

•  Buyer  would  like  to  operate  in  a  decentralized  supply  chain 

with  V  =  0  and  set:  (Q*d,  R*d)  =  TCb(Q,R,  0) 

■  *  Solution:  =  455,  R*d  =  1019,  TC f  =  TCb(Q*d,R*d,  0)  = 

4406.4 

•  Given  buyer’s  (Qd,Rd),  supplier  finds  Nd  —  1  to  minimize  its 
costs. 


i — *  TCj  =  TCs(Q*d,  N*[,  0)  =  3298 

h-  TCd  =  TC$  +  TCd  =  4406.4  +  3298  =  7704.4 
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Finding  Optimal  Policies 


•  Supplier  would  like  to  operate  in  a  centralized  supply  chain 
i  *  Define:  TCC(Q,R,N )  =  TCb(Q,  R,  0)  +  TCS(Q,  N,  0) 

(Q*,  R*c,  N*)  =  TCC  (Q,  R,  N) 

i — *  Solution:  Q*  =  909,  R*  =  1004,  N*  =  1 

i  *  TC£  =  4955.9,  TC'f  =  1649.9,  TCC  =  6605.8 

i — *  Coordination  can  save:  TC+  =  TCd—TCc  —  7704.4-6605.8 
1098.5 
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Observations 


•  Buyer  prefers  a  decentralized  supply  chain:  TCfj  <  TC% 

•  Supplier  prefers  a  centralized  supply  chain:  TC§  >  TCCS 

i — *  Centralized  policy  requires  buyer  to  raise  order  quantity  by 
Q*  -  Q*  =  454  and  can  save  TC+  =  1098.5. 

i — ^  Centralized  no.  of  orders:  Y/Q *  «  11 

i — ^  Seller  offers  buyer  rebate  each  cycle: 


V  =  0.5  • 


TC+ 

y/oi 


1098.5 

0.5 - 

11 


50 


Parties  can  agree  on  other  split  of  TC+ 
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Alternate  Solution  (Chaharsooghi  and  Heydari 

(2010)  —  CH) 

•  Approximate  the  density  function  fx  for  X  using  a  normal 
distribution  using  sample  means  and  variances. 

•  For  this  example,  E(X )  =  fix  —  417.2  and 

Var(X)  =  a\  =  19881 


•  Define:  k  =  (R  —  nx)/crx,  so  that 

(z  —  k)  •  dz  . 

v  27T 

•  CH  solve  analytically  for  a  partial  solution  for  (Q,k,N)  stated 
in  terms  of  the  cost  parameters  and  normal  CDF. 


roc 

S(R)  =  Sk(k)  -ax  —  a X  '  Jk 
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Alternatives  —  MOP  and  CH  Solution 


i — >  We  should  measure  the  effectiveness  of  the  models  by  their 
“value  in  use" 


•  Consider  two  models:  1)  MOP  model;  and  2)  CH  model. 
We  will  compare  the  solutions  obtained  from  the  two  models  by 
simulating  from  the  (unknown)  underlying  “actual"  normal  daily 
demand  and  discrete  LT  distributions  (the  “actual"  model). 
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Evaluating  the  Models 


Decentralized 

Q* 

R* 

N* 

Buyer 

TC 

%  Dec. 

CPU 

(sec) 

MOP  (mixture  dist.) 

455 

1019 

1 

4455 

2.0% 

1.77 

CH  (normal) 

558 

999 

1 

4531 

0% 

0.06 

sc 

CPU 

Coordinated 

Q* 

R* 

N* 

TC 

%  Dec. 

(sec) 

MOP  (mixture  dist.) 

909 

1004 

1 

6628 

3.7% 

2.16 

CH  (normal) 

1012 

926 

1 

6872 

0% 

0.09 

•  The  costs  are  directly  comparable  —  calculated  by  inserting  the 
MOP  and  CH  solutions  into  the  simulation  model  and  running 
100,000  trials. 
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Air  Force  Example  —  F-15/16  Power  Supply 


•  1827  observations  of  daily  demand  (2008-2012)  with  d  =  0.63 
and  —  1.09. 

i — ^  Mode  =  0  (1158  observations);  Maximum  =  8. 

•  100  randomly  sampled  requisitions:  I  =  10.4  and  sj  —  162.8 
(Min  =  l;  Max=73). 

•  Annual  unit  holding  cost:  15%;  Unit  Price:  $224,392; 

—  $33,658;  K ^  =  $5.20. 

•  Annual  unit  shortage  cost  (7r)  -  if  one  unit  short  one  officer  at 
captain  pay  is  50%  productive,  7r=25000. 
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Air  Force  Example 


•  LTD  Distribution  —  MOP  model  and  a  normal  approximation 


•  Use  this  distribution  to  find  optimal  Q*f  R *,  and  N*  policies. 


30 


Air  Force  Example 


•  Comparison  of  solutions  with  MOP  and  normal  approximations 


Holding 

Order 

Shortage 

Total 

Q* 

R* 

N* 

Cost 

Cost 

Cost 

Cost 

Decentralized 

MOP 

13 

15 

1 

2,642,476 

463 

6,438 

2,649,377 

Coordinated 

MOP 

14 

15 

1 

432 

ffUSSl 

Decentralized 

Normal 

7 

22 

1 

858 

274 

mnvzmnm 

Coordinated 

Normal 

9 

21 

1 

3,242,049 

666 

411 

3,243,126 

•  Buyer  costs  calculated  by  implementing  the  policy  with  the 
actual  demand  data  for  2008-2012  (before  considering  coordi¬ 
nation  incentives)  and  lead  times  drawn  randomly  for  each  order 
from  the  empirical  distribution. 
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Order  Quantity/Reorder  Point 


Sensitivity  to  Shortage  Cost  Parameter 
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Conclusions 


•  Mixture  distributions  can  be  used  to  model  the  distribution  for 
demand  during  lead  time  using  strictly  empirical  data  with  no 
limits  on  the  underlying  distribution. 

•  By  using  MOP  distributions  estimated  from  B-spline  functions, 
we  can  perform  integrations  required  to  determine  optimal  order 
quantities,  reorder  points,  and  service  levels  in  closed-form. 

•  Next  steps:  creating  models  under  different  sets  of  assump¬ 
tions,  e.g.  a  vendor-managed  inventory  model;  improving  the 
efficiency  of  the  solution  algorithm. 
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